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Abstract 

Fruit ripening is a complex, genetically programmed process that occurs in conjunction with the differentiation of 
chloroplasts into chromoplasts and involves changes to the organoleptic properties of the fruit. In this study, an 
integrative analysis of the transcriptome and proteome was performed to identify important regulators and pathways 
involved in fruit ripening in a spontaneous late-ripening mutant ('Fengwan' orange, Citrus sinensis L. Osbeck) and its 
wild type ('Fengjie 72-1'). At the transcript level, 628 genes showed a 2-fold or more expression difference between 
the mutant and wild type as detected by an RNA sequencing approach. At the protein level, 130 proteins differed by 
1.5-fold or more in their relative abundance, as indicated by iTRAQ (isobaric tags for relative and absolute quantita- 
tion) analysis. A comparison of the transcriptome and proteome data revealed some aspects of the regulation of 
metabolism during orange fruit ripening. First, a large number of differential genes were found to belong to the plant 
hormone pathways and cell-wall-related metabolism. Secondly, we noted a correlation between ripening-associated 
transcripts and sugar metabolites, which suggests the importance of these metabolic pathways during fruit ripening. 
Thirdly, a number of genes showed inconsistency between the transcript and protein level, which is indicative of post- 
transcriptional events. These results reveal multiple ripening-associated events during citrus ripening and provide 
new insights into the molecular mechanisms underlying citrus ripening regulatory networks. 
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Introduction 

Citrus fruits are among the most important and widely grown 
commodity fruit crops and have a non-climacteric fruit matu- 
ration character (Paul et al, 2012). The growth and develop- 
ment of citrus fruits can be split into three stages: cell division, 
expansion, and ripening (Bain, 1958). At the ripening stage, 
carotenoids, sugars, and other soluble solids are accumulated, 



chlorophyll and organic acid contents are reduced, the cell 
wall is extensively modified, and the concentration of a num- 
ber of volatiles increases (Katz et al. , 201 1; Yu et al. , 2012). 

Plant hormones are very important for fruit development 
and ripening (Theologis, 1992; Fan et al, 1998; Soto et al, 
2013). Recent studies of strawberry (Fragaria ananassa) 



Abbreviations: ABA, abscisic acid; COG, Clusters of Orthologous Groups; DAF, days after flowering; DEG, differentially expressed genes; FDR, false discovery rate; 
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showed that abscisic acid (ABA) promoted fruit ripening and 
played an important role in the regulation of fruit develop- 
ment and ripening (Ji et al, 2012; Jia et al, 2011, 2013b). Jia 
et al. (2011) applied the VIGS (Tobacco rattle virus-induced 
gene silencing) technique to silence a 9-c«-epoxycarotenoid 
dioxygenase gene (FaNCEDl) in strawberry fruit, which is 
key to ABA biosynthesis, resulting in a significant decrease 
in ABA levels and uncoloured fruits. Interestingly, a similar 
uncoloured phenotype was observed in the transgenic RNA 
interference (RNAi) fruits, in which the expression of a puta- 
tive ABA receptor gene encoding the magnesium chelatase 
H subunit (FaCHLHIABAR) was downregulated by virus- 
induced gene silencing and, more importantly, downregula- 
tion of the FaCHLHIABAR gene in the RNAi fruit altered 
both ABA levels and sugar content as well as a set of ABA- 
and/or sugar-responsive genes. Koyama et al. (2010) found 
that ABA application rapidly induced the accumulation of 
anthocyanin and flavonol in berry skins of the Cabernet 
Sauvignon grape, which proved that ABA could stimulated 
ripening and ripening-related gene expression. Recently, 
ABA level was investigated in cucumber ( Cucumis sativus L.), 
and a peak in ABA level was found in pulp before the fruit 
get fully ripe; real-time PCR analysis revealed that the ABA 
content may be regulated by its biosynthesis (CsNCEDs), 
catabolism (CsCYP707Al), and reactivation (CsBGs) genes 
(Wang et al, 2013). Similarly, a study on a spontaneous 
fruit-specific ABA-deficient mutant sweet orange, 'Pinalate', 
showed that ABA was an important regulator for the onset of 
fruit degreening and carotenoid biosynthesis (Romero et al. , 

2012) , although other maturation processes were not affected 
in 'Pinalate' fruit (Rodrigo et al, 2003). Rodrigo et al. (2006) 
found that the content of ABA was increased with the fruit- 
ripening process and that ethylene could stimulate ABA 
content and promote fruit coloration in orange. Recent stud- 
ies have indicated that ABA appears to modulate ripening 
through interference with ethylene and auxin-related genes 
(Zhang et al, 2009; Ji et al, 2012; Sun et al, 2012; Soto et al, 

2013) . ABA was shown to have the ability to control the gene 
expression of the entire system of cell-wall catabolism genes, 
such as XET, PE, and PG (Sun et al, 2012). These studies 
suggest that ABA plays an important role in the regulation 
of fruit ripening. Although ethylene plays a major role in the 
ripening process of climacteric fruits (Theologis, 1992), stud- 
ies on citrus, a non-climacteric group of fruits, have shown 
that ripening-related colour changes in the flavedo portion of 
the fruit peel are regulated by endogenous as well as exog- 
enous levels of ethylene (Goldschmidt, 1997), and the fruit's 
sensitivity to ethylene is also an important factor for fruit 
ripening; even a small amount of ethylene produced by the 
fruit might be sufficient to trigger ripening-related physiologi- 
cal responses in some non-climacteric fruits (Perkins- Veazie 
et al, 1996; Trainotti et al, 2005). 

Sugars were traditionally regarded as the metabolic 
resources required for construction of the carbon skeleton 
and energy supply in plants. Recently, however, many stud- 
ies have suggested that sugars may serve as important sig- 
nals that modulate a wide range of processes in the plant life 
cycle. Notably, intimate connections between sugar and ABA 



signalling have been revealed by the isolation and characteri- 
zation of sugar-insensitive and sugar-hypersensitive mutants 
in Arabidopsis (Smeekens, 2000; Rolland et al, 2006). Other 
studies have demonstrated that ABA and sugars often have 
similar or antagonistic effects on diverse developmental pro- 
cesses (Rolland et al, 2006; Jia et al, 2011). Recent work has 
demonstrated that sucrose is an important signal in the regu- 
lation of strawberry fruit ripening (Jia et al, 2013a). 

The ripening of citrus fruit is accompanied by the synthesis 
of a large number of proteins and the transcription of many 
genes. Several proteins and regulation factors are involved 
in this process. Transcriptomic and proteomic analyses are 
extremely efficient methods for identifying differential expres- 
sion genes at the whole-genome level. Recently, the fruit rip- 
ening of the tomato was investigated at the whole-genome 
level by comparative transcriptomics and proteomics, which 
highlighted the need for combined transcriptomic and pro- 
teomic analyses (Osorio et al, 2011). More specifically, the 
recently developed isobaric tags for relative and absolute 
quantitation (iTRAQ) technology was validated to be power- 
ful in protein profiling (Elias and Gygi, 2007; Gan et al, 2007; 
Noirel et al, 2011; Unwin et al, 2010). However, previous 
studies revealed that the transcripts did not, in fact, always 
coincide well with their final products, the proteins (Osorio 
et al. , 201 1; Pan et al. , 2012). This divergence might be due to 
post-transcriptional and translational processing that regu- 
lates the location, quantity, and efficiency of the proteins in 
the cell. 

In citrus, artificially generated mutants are very difficult 
to generate, but spontaneous mutants are widely found in 
nature. A spontaneous late-ripening mutant from the 'Fengjie 
72-1' orange (Citrus sinensis L. Osbeck), named 'Fengwan', 
has been biochemically characterized (Liu et al, 2006b). 
During natural ripening, the onset of fruit degreening and 
the time of full maturation were delayed in 'Fengwan', com- 
pared with its wild-type (WT) cultivar (Liu et al, 2006b). As 
a result, these two materials, 'Fengjie 72-1' and 'Fengwan', 
have provided a promising platform to investigate the regu- 
lation network of citrus fruit development and ripening. In 
the present study, we examined 'Fengjie 72-1' and 'Fengwan' 
along the ripening periods at the transcriptomic and prot- 
eomic levels. Transcriptional analysis was conducted using 
RNA sequencing (RNA-seq) (Mortazavi et al, 2008), and 
proteomic data were obtained using iTRAQ (Elias and Gygi, 
2007; Gan et al. , 2007; Unwin et al. ,2010; Noirel et al. , 201 1). 



Materials and methods 

Plant materials and RNA preparation 

The WT 'Fengjie 72-1' orange (C. sinensis L. Osbeck) and its spon- 
taneous late-ripening mutant (MT) 'Fengwan' were both cultivated 
in the same orchard (Fengjie, Chongqing City, China). Fruit sam- 
ples were harvested at 150, 170, 190, 210, and 240 d after flowering 
(DAF) from three different trees in 201 1 . Twelve representative fruits 
were sampled from each tree at each developmental stage. After 
separating the pulp from the peel, the pulp and peel were sliced. 
Combining the sliced WT pulp samples with one another (as for the 
MT samples), the samples were frozen rapidly in liquid nitrogen and 
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kept at -80 °C (Liu et al, 2009; Xu et al, 2009). A portion of the 
samples was used for extracting the RNA, as described previously 
(Liu et al, 2006a). Another aliquot was used to extract the protein. 
A portion of the sample powder was used for the determination 
of ABA, sugar, and organic acid composition and concentration. 
Similarly to the sliced pulp samples, the peel samples were frozen 
rapidly in liquid nitrogen and kept at -80 °C for determination of 
the total amount of chlorophyll. 



Analysis of total chlorophyll quantification, soluble sugar, and 
organic acid 

Approximately 0.5 g of peel material was used for determination of 
the total chlorophyll content, according to the method described 
by Li (2000). Each sample was assayed with three replicates. A gas 
chromatography (GC) method was used to determine the soluble 
sugar and organic acid composition and the concentrations of the 
WT and MT fruit pulps harvested at 150, 170, 190, 210, and 240 
DAF, as described previously (Yu et al., 2012). 

Quantification of ABA 

The samples for ABA quantification were prepared according to 
the method described by Pan et al (2010) with some modified. D 6 - 
ABA (Icon Isotopes) was used as an internal standard for ABA. In 
brief, fruit pulp was ground into powder with a mortar and pestle in 
liquid nitrogen, weighed and each sample (200 mg) was transferred 
to 10 ml screw-cap tubes. The samples were kept in liquid nitro- 
gen. Fifty microliters of the working solution of internal standards 
was added to each 10 ml tube containing the frozen plant material 
and 2ml of extraction solvent [2-propanol:H 2 0:concentrated HC1 
(2: 1 :0.002, v/v/v)], was added to each tube. The tubes were placed on 
a shaker at a speed of 200 rpm for 30min at 4 °C. Dichloromethane 



(4 ml) was added to each sample and shaken for 30min at 4 °C. The 
samples were placed into a refrigerated microcentrifuge at 4 °C and 
centrifuged at 13 000g for 5min. Approximately 4 ml of the solvent 
was transferred from the lower phase into a screw-cap vial and the 
solvent mixture was concentrated (not completely dry) using a nitro- 
gen evaporator with nitrogen flow. The samples were redissolved in 
0.2ml of methanol and filtered with 0.22 um organic membrane 
filters. The extracted solution was injected into a reverse-phase 
C 18 Gemini high-performance liquid chromatography (HPLC) 
column for HPLC electrospray ionization tandem mass spec- 
trometry (HPLC-ESI-MS/MS) analysis. An Agilent 1100 HPLC 
(Agilent Technologies), Waters C 18 column (150x2.1 mm, 5 um), 
and API3000 MS-MRM (Applied Biosystems) were used as for the 
analysis. The reaction monitoring conditions (Q1/Q3) of ABA and 
D 6 -ABA were 263.00/153.00 and 269.00/159.10, respectively. Each 
sample was assayed using four replicates. 



RNA-seq, data processing, and gene annotation 

Combining the data of ABA, chlorophyll, soluble sugar, and organic 
acid (Figs 1 and 2), We found that the increase or decrease in sug- 
ars and acids and the largest peel colour change of fruit, as well as 
the ABA peak, all appeared at 170, 190, and 210 DAF, indicating 
that these three ripening stages were the key stages in fruit ripening. 
Thus, the WT and MT fruit pulps harvested at 170, 190, and 210 
DAF were subjected to RNA-seq using an Illumina HiSeq™2000 
at the Beijing Genomics Institute (Shenzhen) in 2011. Briefly, 6 ng 
of the total RNA of each sample was used to enrich the mRNA, 
to construct cDNA libraries, and for sequencing analysis. High- 
quality reads (clean reads) were obtained by removing low-quality 
reads with ambiguous nucleotides, and adaptor sequences were fil- 
tered from the raw reads, and a statistical analysis was conducted to 
summarize the number of clean reads that aligned with the recently 
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Fig. 1 . (A) Phenotypic characterization of fruit development in MT and WT. (B) Chlorophyll content in the fruit peel of MT and WT at five developmental 
stages. (C) ABA content in the fruit pulp of MT and WT at five developmental stages. Bars represent the standard error (SE; n=3). Asterisks represents 
statistically significant differences (*P<0.05; **P<0.01) analysed using Student's t-test (This figure is available in colour at JXB online.) 
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Fig. 2. Trends in organic acid and soluble sugar content in MT and WT fruit pulp (fresh weight) during fruit ripening. Quinic acid (A), citric acid (B), malic 
acid (C), glucose (D), fructose (E), and sucrose (F) content were extracted at 1 50, 1 70, 1 90, 21 0, and 240 DAF. Three individual replicates were used to 
reduce the experimental error. The total organic acid (TOA; G) and total soluble sugar (TSS; H) were the sum of the organic acids (malic acid, citric acid, 
and quinic acid) and the soluble sugars (fructose, glucose, and sucrose), respectively. The ratio of TSS:TOA was determined (I). Bars represent the SE 
(n=3). Asterisks represents statistically significant differences (*P<0.05; **P<0.01) analysed using Student's f-test. 



released reference genome (Xu et al, 2013); gene expression lev- 
els were calculated by the RPKM (reads per kb per million reads) 
method according to Zheng et al. (2012). 

Gene annotation was conducted using the Blast2GO program. 
Gene ontology (GO) enrichment analysis provided all of the GO 
terms that were significantly enriched in differentially expressed 
genes (DEGs) compared with the genome background, and we 
filtered the DEGs that corresponded to biological functions. The 
biological interpretation of the differential genes was further com- 
pleted by assigning them to metabolic pathways using the Kyoto 
Encyclopaedia of Genes and Genomes (KEGG) annotation. 

Sequence alignment and real-time quantitative PCR verification 

Sequence similarity comparisons between Arabidopsis thaliana and 
C. sinensis proteins were performed by BLASTP in the Citrus sinen- 
sis Annotation Project (CAP) (http://citrus.hzau.edu.cn/). A search 
for amino acid sequences of Arabidopsis ABA- and ethylene-related 
proteins was carried out using the National Centre for Biotechnology 
Information (Supplementary Table SI available at JXB online). 

A total of 1 pig of RNA was reverse transcribed for first-strand 
cDNA synthesis using a RevertAid™ First Strand cDNA synthesis 
kit (Fermentas) according to the manufacturer's instructions. The 
gene-specific primer pairs (Supplementary Tables SI and S2 avail- 
able at JXB online), which were designed with the Primer Express 
3.0 software (Applied Biosystems), were used for real-time PCR. 
Reactions were performed with the SYBR Green PCR Master Mix 
in an ABI 7900HT Fast Real-time system. Actin was used as the 
standard to normalize the content of cDNA as described previously 
(Liu et al, 2009). Ten microlitres of the reaction mixture was added 
to each well. The thermal cycling program was set at 50 °C for 2min, 
95 °C for 1 min, and 40 cycles of 95 °C for 15 s and 60 °C for 1 min. 
The output results were analysed by the instrument on-board soft- 
ware Sequence Detector Version 1.3.1 (PE Applied Biosystems). 



The real-time PCR was conducted with five replicates for each sam- 
ple, and data are indicated as means ±standard error (SE) («=3). 

Protein extraction, protein digestion, and iTRAQ labelling 

The WT and MT fruit pulps harvested at 170, 190, and 210 DAF 
were used for protein extraction, which was performed using phenol 
extraction, as described previously (Pan et al., 2009). The protein 
concentration was determined with a Bio-Rad Protein Assay kit 
based on the Bradford method using BSA as a standard. Two inde- 
pendent protein extractions were performed. 

One hundred micrograms of total protein from each sample solu- 
tion was used for protein digestion. First, the protein was digested 
with Trypsin Gold at a protein:trypsin ratio of 20:1 (w/w) at 37 °C 
for 4h. The Trypsin Gold was then added at a proteimtrypsin 
ratio of 20:1 (w/w) once more, and the result was digested for 8h 
continuously. 

After trypsin digestion, the peptide was dried by vacuum centrif- 
ugation. The peptide was reconstituted in 0.5 M triethylammonium 
bicarbonate buffer and processed according to the manufacturer's 
protocol for 8-plex iTRAQ (Applied Biosystems). Briefly, 1 U of 
iTRAQ reagent (defined as the amount of reagent required to 
label 100 iig of protein) was thawed and reconstituted in 70 ul of 
isopropanol. Peptides from digestion were labelled with different 
iTRAQ tags in the same group by incubation at room tempera- 
ture for 2h. Samples of WT fruit pulp harvested at 170, 190, and 
210 DAF were each labelled with iTRAQ reagents with molecular 
masses of 114, 116, and 119Da, respectively. Samples of MT fruit 
pulp harvested at 170, 190, and 210 DAF were each labelled with 
iTRAQ reagents with molecular masses of 113, 115, and 118 Da, 
respectively. After labelling, the peptide mixtures were pooled and 
dried by vacuum centrifugation. The pooled mixtures of iTRAQ- 
labelled peptides were fractionated by strong cationic exchange 
(SCX) chromatography. 
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Fractionation by SCX 

For SCX chromatography using a Shimadzu LC-20AB HPLC 
Pump system, the iTRAQ-labelled peptides were reconstituted with 
4ml of buffer A (25 mM NaH 2 P0 4 in 25% acetonitrile, pH 2.7) and 
loaded onto a 4.6x250 mm Ultremex SCX column that contained 
5 um particles (Phenomenex). The peptides were eluted at a flow rate 
of 1 ml mkr 1 with a gradient of buffer A for lOmin, 5-35% buffer B 
(25 mM NaH 2 P0 4 , 1 M KC1 in 25% acetonitrile, pH 2.7) for 1 1 min 
and 35-80% buffer B for 1 min. The system was then maintained in 
80% buffer B for 3 min before equilibrating with buffer A for 10 min 
prior to the next injection. Elution was monitored by measuring the 
absorbance at 214nm, and the fractions were collected every 1 min. 
The eluted peptides were pooled as 12 fractions, desalted using a 
Strata X CI 8 column (Phenomenex), and vacuum dried. 

LC-ESI-MS/MS analysis by LTQ Orbitrap HCD 

Each fraction was resuspended in a certain volume of buffer A (2% 
acetonitrile, 0. 1% formic acid) and centrifuged at 20 OOOg for 10 min. 
In each fraction, the final concentration of peptide was approxi- 
mately 0.5 p.g uT 1 , on average. A total of 10 of supernatant was 
loaded onto a Shimadzu LC-20AD nanoHPLC by the autosampler 
onto a 2cm CI 8 trap column (inner diameter, 200 um), and the 
peptides were eluted onto a resolving 10cm analytical C18 column 
(inner diameter, 75 um) made in house. The samples were loaded 
at 15 ill mkr 1 for 4 min; the 44 min gradient was then run at 400 nl 
min -1 from 2 to 35% B (98% acetonitrile, 0. 1% formic acid), followed 
by a 2 min linear gradient to 80%, maintenance at 80% B for 4 min, 
and finally a return to 2% for 1 min. 

The peptides were subjected to nanoESI followed by MS/MS 
in an LTQ Orbitrap Velos (Thermo Scientific) that was coupled 
online to the HPLC. Intact peptides were detected in the Orbitrap 
at a resolution of 60 000. Peptides were selected for MS/MS using 
a high-energy collision dissociation (HCD) operating mode with a 
normalized collision energy setting of 45%; the ion fragments were 
detected in the LTQ. A data-dependent procedure that alternated 
between one MS scan followed by eight MS/MS scans was applied 
for the eight most abundant precursor ions above a threshold ion 
count of 5000 in the MS survey scan with the following Dynamic 
Exclusion settings: repeat counts, 2; repeat duration, 30 s; and exclu- 
sion duration, 120 s. The electrospray voltage applied was 1.5kV. 
Automatic gain control was used to prevent overfilling of the ion 
trap; 1 x 10 4 ions were accumulated in the ion trap for the genera- 
tion of the HCD spectra. For the MS scans, the mlz scan range was 
350-2000 Da. 



Database search and quantification 

Mascot software version 2.3.02 (Matrix Science) was used to simul- 
taneously identify and quantify the proteins. In this version, only 
unique peptides used for protein quantification can be chosen, which 
provides a more precise quantification of the proteins. Searches were 
made against the database (44272 sequences; http://citrus.hzau.edu. 
cn/orange/) (Xu et al, 2013). Spectra from the 12 fractions were 
combined into one MGF (Mascot generic format) file after loading 
the raw data, and the MGF file was searched. The search parameters 
were as follows: trypsin was chosen as the enzyme with one missed 
cleavage allowed; fixed modifications of carbamidomethylation at 
Cys, variable modifications of oxidation at Met; peptide tolerance 
was set at 0.05 Da, and MS/MS tolerance was set at 0.1 Da. The pep- 
tide charge was set to 'Mr', and the monoisotopic mass was chosen. 
An automatic decoy database search strategy was also employed to 
estimate the false discovery rate (FDR); the FDR was calculated as 
the number of false-positive matches divided by the total number 
of matches. In the final search results, the FDR was less than 1.5%. 
The iTRAQ 8-plex was chosen for quantification during the search. 

The search results were passed through additional filters before 
exporting the data. For protein identification, the filters were set as 
follows: significance threshold P<0.05 (with 95% confidence) and an 



ion score or expected cut-off of less than 0.05 (with 95% confidence). 
For protein quantitation, the filters were set as follows: 'median' was 
chosen for the protein ratio type (http://www.matrixscience.com/ 
help/quant_config_help.html); the minimum precursor charge was 
set to 2+, and the minimum peptide was set to 2; only unique pep- 
tides were used to quantify the proteins. The median intensities were 
set to normalization, and outliers were removed automatically. The 
peptide threshold was set as above for identity. 

Results 

Changes in soluble sugars, organic acid, and 
chlorophyll content during fruit ripening in MTand WT 

Fruit taste is determined by the content and types of soluble 
sugars and organic acids. Thus, the levels of soluble sugars 
and organic acids of MT and WT were determined during 
fruit ripening. In ripe citrus fruit, the main soluble sugars are 
sucrose, fructose, and glucose. The trends for organic acid and 
soluble sugar content are shown in Fig. 2. The content of the 
soluble sugars increased markedly during the stages of fruit 
ripening in both WT and MT (Fig. 2D-F, H). The content of 
the soluble sugars in MT were all obviously lower than that 
in WT from 170 to 240 DAF (Fig. 2D-F). In contrast, the 
concentrations of citric acid, quinic acid, and total organic 
acid (TOA) decreased with the development of fruit ripening 
in both WT and MT (Fig. 2A, B, G). However, the malic acid 
concentrations increased during the stages of fruit ripening 
(Fig. 2C). In addition, we observed that the concentrations 
of citric acid, malic acid, and TOA were higher in MT com- 
pared with WT from 170 to 240 DAF, but the concentrations 
of quinic acid were lower. Thus, the decrease in TOA, citric 
acid, and quinic acid in MT and WT fruit and the increase 
in total soluble sugar (TSS), malic acid, and sucrose, and the 
change in the ratio of TSS:TOA were significant after 170 
DAF. The chlorophyll content of the WT and MT peels were 
also determined at five stages (Fig. IB). We observed that, at 
stages following mature green (170 DAF), the mutant fruit 
diverged in its development from the WT, as the mutant is 
late ripening (Fig. lA).This finding suggested that 170 DAF 
was the turning point at which the MT fruit diverged in its 
development from the WT. 

Fruit transcriptome difference between MT and WT 
during fruit ripening 

One of the primary goals of transcriptome sequencing 
(Supplementary Results available at JXB online) is to com- 
pare the gene expression levels in two genotypes. In this study, 
we used a stringent value of FDR <0.001 and a P value <0.05 
as the threshold to judge the significant differences in the gene 
expressions. A total of 628 genes were differently expressed 
(>2-fold) between MT and WT (Supplementary Table S3 
available at JXB online); of these, 352 were at 170 DAF, 270 
were at 190 DAF, and 90 were at 210 DAF (Supplementary 
Table S4 available at JXB online). At all three developmen- 
tal stages, the number of upregulated genes was less than the 
number of downregulated genes, and the number of DEGs 
was decreased following fruit development (Fig. 3). 
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Fig. 3. Number of DEGs between MT and WT at 1 70, 1 90, and 21 0 DAF. 

To gain insight into the functional categories that were 
altered between MT and WT, GO categories were assigned 
to the 628 DEGs with a BLASTX hit using Blast2GO. Fig. 4 
shows the distributions of the GO terms according to the 
corresponding biological process, molecular function, and 
cellular component. Cell (99 genes distributed in 170 DAF, 
78 genes distributed in 190 DAF, and 24 genes distributed in 
210 DAF), intracellular (60, 32, and 11), organelle (57, 31, 
and 1 1) and membrane (39, 29, and 8) were the major catego- 
ries annotated for the cellular component (Fig. 4A). Binding 
(86, 60, and 18), catalytic activity (118, 88, and 27), hydrolase 
activity (36, 36, and 8), transferase activity (43, 25, and 7), 
and oxidoreductase activity (36, 28, and 10) were the major 
categories annotated under molecular function (Fig. 4B). 
According to the GO terms of the biological process, the 
majority of DEGs appeared to be related to four major bio- 
logical changes, including the primary metabolic process (68, 
39, and 12), cellular process (71, 41, and 16), cellular meta- 
bolic process (62, 31, and 12), and response to stimulus (41, 
18, and 9) (Fig. 4C). In addition, we noted that the number of 
DEGs distributed to the secondary metabolic processes (7, 6, 
and 4) was much less than the number distributed to the pri- 
mary metabolic processes (68, 39, and 12) (Fig. 4C). KEGG 
analysis assigned the differential genes to 54 metabolic path- 
ways (each of which contained three or more differential 
genes). The complete list of metabolic pathways is provided 
in Supplementary Table S5 available at JXB online. Table 1 
lists the metabolic/biological pathways that contained over 
five differential genes. Notably, 16 differential genes were pre- 
dicted to be involved in plant hormone signal transduction at 
170 DAF. Similar to previous studies (Zhang et al. , 2009; Sun 
etal, 2012; Jia et al, 2013a), the metabolic pathways of plant 
hormone signal transduction, carotenoid biosynthesis, flavo- 
noid biosynthesis, and starch and sucrose metabolism were 
involved in fruit ripening (Table 1). 

Plant hormones are very important for fruit development 
and ripening. In this study, as shown in Table 2, 11 genes 
were found to be involved in plant hormone metabolism, 
including five genes involved in ABA synthesis and signal 
transduction (Cs5gl4370, orangel.lt03687, Cslg09250, 
Cs6gl9380, and Cs8gl3770) and six genes related to ethylene 
(Cs2g02500, Cs4gl3870, Cs5g29870, Cs9g08850, Cs4gl7960, 
and Cs8gl2880), and most of these genes were downregu- 
lated in MT. Notably, some key enzymes for synthesizing 
plant hormones were found among these 1 1 genes, such as 



NCED1 (Cs5gl4370), ACO (Cs2g02500, Cs4gl3870), AAO 
(Cs8gl3770), and ABAS'-hydroxylase (orangel.lt03687, 
Cslg09250, and Cs6gl9380) (Table 2). Some important regu- 
lators involved in ethylene signal transduction were identified, 
including ETR (Cs9g08850, Cs8gl2880), EBF4 (Cs4gl7960), 
and ERF1B (Cs5g29870) (Table 2). In addition, 15 genes 
were identified to be involved in cell-wall metabolism, includ- 
ing five XETs, five PEs, two PGs, and three other genes 
(Table 2). The five XETs displayed a mixed expression pat- 
tern (two genes upregulated and three genes downregulated), 
and (3-glucosidase also displayed a mixed expression pattern 
(Table 2). However, the five PEs were all downregulated in 
MT. In contrast, the two PGs were all upregulated in MT. The 
three sucrose-related genes were identified as sucrose synthase 
(Cs6gl5930), sucrose-phosphate synthase (Cs5gl9060), and 
citrus sucrose transporter 1 (Cs3g22560); these genes were all 
downregulated in MT (Table 2). 

To obtain a global overview of the genes involved in fruit 
ripening, we summarized the change in the gene expression 
levels during fruit ripening for WT and MT (Supplementary 
Tables S7 and S8 available at JXB online). We observed that 
there were 1036 developmentally DEGs in WT and 1406 
in MT. KEGG analysis assigned the DEGs of WT and 
MT to different metabolic pathways. Supplementary Table 
S6 available at JXB online lists some important metabolic/ 
biologic pathways. Notably, most of the DEGs of WT and 
MT were annotated to 'Plant-pathogen interaction', 'Plant 
hormone signal transduction', 'Stilbenoid, diarylhepta- 
noid and gingerol biosynthesis', 'Microbial metabolism in 
diverse environments', and 'Phenylpropanoid biosynthesis' 
(Supplementary Table S6 available at JXB online). 

Some of the important DEGs are listed in Table 3. 
Several DEGs were associated with plant hormone bio- 
synthesis and signal transduction (Table 3). A gene that 
encoded NCED1 (Cs5gl4370) and two PP2Cs (Cs8g20420, 
Cs9g 18020) showed an increase during fruit ripening. An 
NCED2 (Cs8gl4150), an SAMD (Cs4g02260), and two 
ABA S'-hydroxylase (Cs8g05940, Cs3g23530) genes showed 
a decrease during fruit ripening. Some of the genes that 
are involved in ABA or ethylene signal transduction dis- 
played a different expression pattern; for example, ERF1B 
(Cs5g29870) and PYL4 (Cs7g30500) showed a decrease 
during fruit ripening; however, PP2Cs showed an increase 
during fruit ripening (Table 3). Some of the genes involved 
in sucrose and cell-wall metabolism also underwent a large 
change during fruit ripening, including those coding for 
sucrose synthase (Cs6gl5930), sucrose-phosphate synthase 
(Cs5gl9060), citrus sucrose transporter 1 (Cs3g22560), five 
XETs (Cs4g03200, Cs4g03220, Cs4g03130, Cs4g03060, 
Cs4g03210), seven pectinesterase genes (orangel.lt02719, 
Cs4gl5560, Cs4g06650, Cs4g06630, Cs3gl4650, Cs3gl4610, 
Cs2g07660) and two PGs (orange l.lt00402, Cs9g03730), 
among others (Table 3). Some homologues of (3-glucosidase 
were identified in this study; these homologues displayed 
different expression patterns (Tables 2 and 4). In addition, 
four a-l,4-galacturonosyltransferase (GAUT) genes were 
differentially expressed among different ripening stages in 
both WT and MT and exhibited decreased expression during 
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Fig. 4. Functional categorization of genes with significant transcriptional differences between MT and WT. The genes were categorized based on GO 
annotation, and the number of each category is displayed based on cellular components (A), molecular function (B), or biological process (C). 



fruit ripening (Table 3). Interestingly, 1 1 heat-shock proteins 
(HSPs) were identified, most of which were upregulated 
at 190 DAF in MT; however, these HSPs underwent little 
change in WT (Table 3). Another four L-ascorbate oxidase 
genes (Cs2g21220, Cs5g09230, Cslg23090, and Cs2g29090) 
were identified, which were all decreased during fruit ripening 
in both WT and MT (Table 3). 



Verification of differentially expressed genes during fruit 
ripening 

Transcriptional regulation revealed by RNA-seq data was 
confirmed in a biologically independent experiment using 
a real-time quantitative reverse transcription PCR. A total 
of 19 genes were selected to design gene-specific primers 
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Table 1. KEGG pathways of more than five differential genes in MT compared with WT at 170, 190, and 210 DAF 



KEGG pathway 


Gene number 








4 -7f\ nAC 
1 /U UAr 


iyu UAr 


£.i\j UAr 


Plant hormone signal transduction 


16 


1 1 


1 


Phenylpropanoid biosynthesis 


15 


13 


5 


Microbial metabolism in diverse environments 


15 


9 


4 


Stilbenoid, diarylheptanoid, and gingerol biosynthesis 


14 


8 


4 


Plant-pathogen interaction 


13 


14 


2 


Protein processing in endoplasmic reticulum 


12 


1 


- 


Flavonoid biosynthesis 


11 


7 


4 


Polycyclic aromatic hydrocarbon degradation 


10 


6 


3 


Starch and sucrose metabolism 


10 


7 


2 


Limonene and pinene degradation 


9 


6 


3 


Aminobenzoate degradation 


9 


7 


3 


Pentose and glucuronate interconversions 


8 


3 


1 


Carotenoid biosynthesis 


6 


- 


2 


Glycine, serine and threonine metabolism 


6 


2 


2 


Phenylalanine metabolism 


6 


8 


1 


Antigen processing and presentation 


5 


1 


- 


Cyanoamino acid metabolism 


5 


2 


- 


Apoptosis 


5 


3 




Diterpenoid biosynthesis 


4 


5 


1 


Methane metabolism 


4 


5 


2 


Lysosome 


3 


5 


3 


Glycerophospholipid metabolism 


3 


5 


3 


RNA transport 


1 


5 


1 



(Supplementary Table S2 available at JXB online) for real- 
time PCR analysis (Fig. 5). A linear regression analysis 
showed an overall correlation coefficient of R-0.193, which 
indicates a good correlation between transcript abundance 
assayed by real-time PCR and the transcription profile 
revealed by RNA-seq data (Fig. 5C). 

Transcriptional regulation of ABA, ethylene 
biosynthesis, and signal transduction related genes 
during orange fruit development and ripening: influence 
of endogenous ABA levels 

Through the analysis of the RNA-seq data, we found that 
hormone-related genes in the orange fruit ripening changed 
greatly, especially the genes of ABA biosynthesis and signal 
transduction pathways, followed by the genes of ethylene 
biosynthesis and signal transduction pathways. In previous 
studies, it has been shown that ABA and ethylene in the mat- 
uration process play important roles, and between them there 
is a very close interaction. To understand better the transcrip- 
tional regulation of ABA and ethylene-related genes in fruit 
ripening, and their relationship with the endogenous ABA 
levels, expression analysis of the 18 ABA-related genes and 
1 5 ethylene-related genes was carried out together with meas- 
urement of ABA in the pulp of fruits of WT and MT. 

Fruits harvested at 150, 170, 190, 210, and 240 DAF were 
selected. As expected, the pulp of WT fruits reached the high- 
est ABA levels at 170 DAF, while the MT fruits reached the 
highest ABA levels at 190 DAF, 20 d later than WT. Although 
it took WT and MT a different number of days to reach the 



highest level of ABA, the final concentration of the highest 
ABA levels of these cultivars was almost equal (Fig. 1C). In 
the pulp of full-colour fruits from both cultivars, an impor- 
tant decrease in ABA content was observed, but levels in the 
WT fruit remained higher than that in the MT (Fig. 1A, C). 

For the differential ABA accumulation in WT and MT 
during ripening, the analysis of ABA-related genes revealed 
a differential regulation between both cultivars, but the 
overall expression patterns of most of genes were the same, 
and the expression of some genes in MT peaked 20 d later 
than that in WT, such as ABA 8' '-hydroxylase 1, CsHAI, 
and CsHABl (Fig. 6A, B). NCED is a key gene in the bio- 
synthesis of ABA, we analysed CsNCED2 and CsNCED4, 
and found that the expression of these two genes in WT 
were higher than in MT at the stage of 150 DAF but had 
no differences at the other four stages. The expression of 
CsNCEDl in WT was higher than that in MT through- 
out the five stages, and the expression of this gene took 
on a continuously upregulated trend in the two cultivars. 
CsAAO, ABA 8' '-hydroxylase 1, and ABA 8' '-hydroxylase 3 
transcript levels fluctuated during ripening in WT and MT 
and had significant differences between these two cultivars, 
especially the ABA8 '-hydroxylase 1 gene, which showed a 
delayed expression in MT (Fig. 6A). The expression patterns 
of CsPYL2, CsPYL8, CsPYL9, CsABI2, and CsSNRK2.2 
in WT and MT were the same, but their transcript levels 
showed differences (Fig. 6B); CsPYRl, CsHAIl, CsAHG3, 
and CsHABl differed in expression patterns and expres- 
sion levels between WT and MT, and the expression of 
these genes in MT showed a typical delay, which was very 



Ripening regulators identified through transcriptome and proteome | 1659 

Table 2. A list of some important genes that are differentially expressed between MTand WT, and are involved in plant hormone 
biosynthesis/signal transduction, cell-wall metabolism and sucrose-related genes 

The BLAST annotation includes the gi number, e-value, gene description and species for the gene. 



Gene ID 


log 2 ratio (MT/WT) 




BLAST annotation 




170 DAF 


190 DAF 


210 DAF 




Plant hormone biosynthesis and signal transduction 




L/ojy 1 M-O / u 


-2 21 






gi 68300870/0/9-cis-epoxycarotenoid dioxygenase 1 [Citrus sinensis] 


ur d.1 iyc I . I lUuuu I 


-2.12 






gi 255544242/1 .64691 e-52/ABA 8'-hydroxylase [Ricinus communis] 


i yuyzou 


I .O^l 






gi 225445688/1 .75384e-79/ ABA 8'-hydroxylase [Vitis vinifera] 


L-Sogi yooU 


— 1 .U/ 






gi|367465454/0/ABA 8'-hydroxylase [Citrus sinensis] 


Os8g13770 






-1 .08 


gi 147841 197/0/abscisic-aldehyde oxidase [Vitis vinifera] 


os^gu^oUU 


-1.18 






gi 642062/2. 15051 e-87/1 -aminocyclopropane-1-carboxylate oxidase [Pelargoniumxhortorum] 


Os4g 13870 




-1 .74 


-2.27 


gi|359484312/3.38653e-128/PREDICTED: 1 -aminocyclopropane-1 -carboxylate oxidase 1 [Vitis vinifera] 


Os5g29870 


-2.47 






gi 224082796/2. 90568e-45/ Ethylene-responsive transcription factor 1 B [Populus trichocarpa] 


Gs9g08850 


-2.12 






gi|283520944/0/ethylene response 2 [Citrus sinensis] 


Gs4g1 /960 


-1 .94 






gi|224090409/0/ein3-binding f-box protein 4 [Populus trichocarpa] 


Cs8g 12880 


-1.13 






gi|283520946/0/ethylene response 3 [Citrus sinensis] 


Cell-wall metabolism 










OS2g1 4920 


-1 .93 






gi|1008904/9.98188e-140/xyloglucan endotransglycosylase [Tropaeolum majus] 


Gs8g03o50 


-1 .62 






gi|356534254/1 . 1 001 5e-1 33/xyloglucan endotransglucosylase/hydrolase protein 1 0-like [Glycine max] 


Gs8g1 5/20 


-1 .14 


-2.52 




gi|156739650/3.09167e-61/xyloglucan endotraglucosylase/hydrolase [Vigna angularis] 


Cs4g03220 




1 .57 




gi|225446121/3.36075e-123/xyloglucan endotransglucosylase/hydrolase protein 23 [Vitis vinifera] 


Cs1g21130 




1 .29 




gi|308229788/3.3281 3e-94/ xyloglucan endotransglucosylase/hydrolase [Gossypium hirsutum] 


Cs7g19180 


-2.23 






gi|356499962/1 . 1 0841 e-1 05/pectinesterase PPME1 -like [Glycine max] 


Cs4g1 5560 


-1.88 






gi|255546301/8. 86701 e-60/Pectinesterase PPE8B precursor [Ricinus communis] 


Cs9g14330 


-1.76 




-1.13 


gi|255539737/3.46069e-48/pectinesterase [Ricinus communis] 


Cs5g09370 


-1.03 






gi|255537037/9.95994e-65/pectinesterase [Ricinus communis] 


Cs3g1 4650 




-10.48 




gi|2241 23034/7. 78899e-118/pectinesterase [Populus trichocarpa] 


Cs9g 14590 


-2.91 






gi|1 15548295/0/beta-fructofuranosidase [Citrus sinensis] 


Cs9g01630 


1.41 


1.56 




gi|2241 01 497/9.34362e-81 /predicted protein [Populus trichocarpa] 


Cs7g08820 


1.08 


1.03 




gi|255568780/2.94863e-76/polygalacturonase, putative [Ricinus communis] 


Cs8g03370 


1.45 


1.24 




gi|356523324/1.43914e-50/PREDICTED: beta-glucosidase 11 -like [Glycine max] 


Cs8g08680 




-1.16 




gi|255573163/0/Beta-glucosidase, putative [Ricinus communis] 


Sucrose-related genes 








Cs6g 15930 




-2.67 




gi|3571 23064/0/ sucrose synthase 2-like [Brachypodium distachyon] 


Cs5g19060 




-1.97 




gi|255561 468/1 .91 935e-1 67/sucrose-phosphate syntase, putative [Ricinus communis] 


Cs3g22560 




1.42 




gi|21 06392 1/0/citrus sucrose transporter 1 [Citrus sinensis] 



consistent with the late-ripening characteristic of MT. 
The expression of CsPYL4, CsAHGl, and CsABIl in WT 
were 2-fold higher than that in MT at 150 DAF stage, but in 
the other stages (190-240 DAF), these genes showed minor 
differences between both cultivars (Fig. 6B). 

The analysis of the ethylene-related genes revealed a mixed 
expression pattern between both cultivars. Some gene expres- 
sion patterns were consistent between both cultivars, such as 
CsCTRl and CsEIL3, while other gene expression patterns 
were largely different between both cultivars, such as CsASC2 
and CsMPK6 (Fig. 6C). ACOl, ASCI, and ASC2 are the 
key genes for the biosynthesis of ethylene. CsACOl was con- 
tinuously upregulated in expression in both cultivars, and the 
expression levels in WT were higher than that in MT at 190 
and 240 DAF (Fig. 6C). CsASCl and CsASC2 in MT reached 
maximum expression levels at 170 DAF At 150, 210, and 240 
DAF, the transcript levels of these two genes in MT were lower 
than those in WT (Fig. 6C). The expression patterns of minor 
genes of the ethylene signal transduction pathway were simi- 
lar, but the majority of genes had differences between WT and 



MT, such as CsERSl, CsEIN3, CsEIN2, CsEIL4, CsERFl, 
and CsEBFl. CsMPK6 and CsERF2 had differences between 
WT and MT, and the expression level of CsMPK6 in WT was 
13-fold higher than that in MT at 150 DAF The transcript 
level of CsERF2 in MT reached a maximum at 210 DAF stage 
and it was 2-fold higher than that in WT, but the expression of 
CsERF2 was higher in WT at 1 50, 170, and 240 DAF (Fig. 6C). 

Proteomic analysis 

In a parallel analysis, a comparative proteome survey was 
performed on MT and WT by the iTRAQ technique to 
complement the transcriptome study. A total of 1566 and 
1501 proteins were identified in the two biological experi- 
ments, respectively (Supplementary Table S9 available at JXB 
online). For both of the biological repeats, the FDR was less 
than 0.67% (0.47% for experiment 1 and 0.67% for experi- 
ment 2). A total of 1839 proteins were identified by merging 
the data obtained from the two biological replicates, with an 
overlap of more than 78% (Fig. 7A). 
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Table 3. A list of some of the important differentially expressed genes between the different ripening stages in MT and WT, including 
genes involved in plant hormone biosynthesis/signal transduction, cell-wall metabolism, sucrose-related genes, ascorbate and aldarate 
metabolism, and protein processing in the endoplasmic reticulum 

Here, 170, 190 and 210 indicate 170, 190 and 210 DAF, respectively. The change fold is shown as a log 2 ratio. 



Gene accession Gene description WT fold change MT fold change 

190 vs 170 210 vs 190 190 vs 170 210 vs 190 



Plant hormone signal transduction 



Cs8g20420 


PP2C 25 


1.25 


1.37 




2.20 


Cs9g1 8020 


PP2C8 


1.04 




1.46 




Cs4g1 7960 


ein3-binding f-box protein 4 


-1.98 








Cs9g08850 


Ethylene response 2 


-1.96 


1.44 




1.99 


Cs7g30500 


PYL4 


-1.77 




-1.32 




Cs4g1 8640 


PYL2 






11.21 




Cs3g1 8000 


Serine/threonine-protein kinase SRK2 






1.24 




Cs1g03300 


EREBP-like factor 






1.25 




Cs8g 12880 


Ethylene response 3 








1.01 



Plant hormone biosynthesis 

Cs5g14370 

Cs8g05940 

Cs1g09250 

Cs8g14150 

Cs3g23530 

Cs7g 14820 

orange1.1t00416 

Cs9g01410 

Cs4g 13870 

Cs4g02260 

Cs2g20590 

Cs2g02500 
Cell-wall metabolism 

Cs4g03200 

Cs4g03220 

Cs4g03130 

Cs4g03060 

Cs4g03210 

orange1.1t02719 

Cs4g 15560 

Cs4g06650 

Cs4g06630 

Cs3g 14650 

Cs3g14610 

Cs2g07660 

orange1.1t00402 

Cs9g03730 

Cs2g22990 

Cs7g 16980 

Cs5g27650 

Cs3g24380 

Cs9g02570 

Cs8g08680 

Cs7g23860 

Cs7g01340 
Sucrose-related genes 

Cs6g 15930 

Cs5g 19060 

Cs3g22560 



NCED1 

ABA 8'-hydroxylase 
ABA 8'-hydroxylase 
NCED2 

ABA 8'-hydroxylase 

NCED4 

ACS 

SAMD 

AC01 

SAMD 

ACO 

ACO 

XET 
XET 
XET 
XET 
XET 

Pectinesterase 3 
Pectinesterase PPE8B 
PME3 

Thermostable pectinesterase 

Pectinesterase 

Pectinesterase 

Pectinesterase 2.2 

Polygalacturonase 

Polygalacturonase 

GAUT1 

GAUT1 

GAUT1 

GAUT8 

p-Glucosidase 

p-Glucosidase 

p-Glucosidase 

p-Glucosidase 

Sucrose synthase 2 
Sucrose-phosphate synthase 
Citrus sucrose transporter 1 



2.69 
-2.44 
-1 .37 



-1.51 



-1.75 

-2.39 
-1.34 



-1.72 

-1.28 

-1.10 
-1.15 
-1.02 
1.92 
1.45 



Ascorbate and aldarate metabolism 
Cs2g21 220 L-Ascorbate oxidase 

Cs5g09230 L-Ascorbate oxidase 



-1.21 



-2.61 
-1.76 
-1.01 
-3.14 

2.80 
-2.53 
-2.58 
-2.41 



-3.34 
-3.28 



-1.14 
-1.44 
-10.48 

-11.14 

-1.32 
-1.74 



-1.51 
-1.47 



-10.81 
-2.61 

-1 .49 



4.40 
-1 .36 



-1 .40 
-1 .37 



3.14 
1.89 



-1.82 
1.98 
-1 .88 



1.12 
-2.31 



2.75 
-1.09 



-1 .38 



-1.63 



1.28 
-2.15 



-2.04 



-2.28 
-1.11 
2.27 
-2.95 

-5.85 

-1.25 
-4.28 

-11.17 
-3.36 
-1 .41 

-3.19 
-1.41 



-2.79 

-1.60 
-1,17 



-1.22 



-1 .44 
-2.18 
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Table 3. Continued 



Cs1g23090 


L-Ascorbate oxidase 


-2.22 


-1 .72 


-1 .59 


-2.47 


Cs2g29090 


L-Ascorbate oxidase 


-2.84 


-3.40 


-3.34 


-10.30 


Cs5g11560 


UDPglucose 6-dehydrogenase 


-1.36 




-2.05 




Protein processing in endoplasmic reticulum 










orange1.1t03636 


HSP20 family protein 






2.44 




orangG i . i XuoZoo 


norzu Tamny protein 






2.97 




Cs8g 19540 


HSP20 family protein 








1.12 


Cs8g 18360 


HSP18.1A 






1.14 




Cs8g 18020 


HSP18.1A 






1.86 




Cs7g29040 


Heat-shock 70kDa protein 1/8 


1.65 






_ 


Cs6g07320 


HSP20 family protein 






1.45 




Cs5g04230 


HSP20 family protein 




-1 1 .03 






Cs4g05880 


HSP20 family protein 






1.99 




Cs2g24360 


HSP20 family protein 






1.49 




Cs2g24040 


HSP20 family protein 






2.04 





Table 4. Correlation between the expression ratios (MT/WT) 

The relative change in abundance (MTA/VT) is shown as a log 2 value from 1 70, 1 90, and 21 0 DAF orange fruit. 



Protein ID 



Proteins (log 2 fold 
change, MT/WT) 



Transcripts (log 2 fold 
change, MT/WT) 



GenBank BLAST annotation 



E-value 



Genbank ID 



Annotation 



1 70 DAF 

Cs7g04020.1 -0.590 

Cs8g1 8020.1 -0.902 

Cs6g09150.1 -0.683 

orange1.1t03414.1 -1.388 

Cs2g24040.1 -1.168 

Cs8g03550.1 -0.665 

190 DAF 

Cs9g01 630.1 0.614 

orange1.1t00340.1 0.932 
210 DAF 

Cs3g12080.1 -1.019 



1.066 
-1 .522 
-1.421 
-1 .267 
-1 .733 
-1 .620 



1.559 
1.019 

-2.379 



8.00E-96 
4.00E-55 
5.00E-62 
6.00E-42 
2.00E-41 
2.00E-147 



0 

1 .00E-53 



1 .00E-1 1 



gi|255540985 
gi|31 593271 8 
gi|255571441 
gi|255563723 
gi|357497003 
gi|356534254 



gi|255556512 
gi|255541538 

gi|2274917 



40S ribosomal protein S9 
HSP18.1A 

Ferritin 

Early nodulin-like protein 1 
18.2kDa class I heat-shock protein 
Xyloglucan endotransglucosylase/ 
hydrolase protein 10-like 

Polygalacturonase precursor 
Remorin 

Cu/Zn superoxide dismutase 



According to previous studies (Zhu et ah, 2009; Lan et al, 
2011), regulated proteins are based on a 1.2- 1.5 -fold change 
threshold. We used a 1.5-fold cut-off to designate changes 
in abundance as significant for the regulated proteins. Using 
these criteria, a total of 130 proteins were classified as being 
differentially expressed between MT and WT (Supplementary 
Table S10 available at JXB online). In general, we observed 
that, across the three stages analysed, fruits at the 170 DAF 
stage showed the highest number of differentially abundant 
proteins in MT compared with WT. This total could be sub- 
divided into 70, 45, and 37 proteins that vary in abundance 
at 170, 190 and 210 DAF, respectively (Fig. 7B). The number 
of downregulated proteins was greater than the number of 
upregulated proteins at the three ripening stages, which was 
consistent with the RNA-seq data (Fig. 3). To obtain func- 
tional information about the 130 differential proteins iden- 
tified, the biological processes and cellular components and 
molecular functions based on the Blast2GO program was 
searched. The results of the biological process categories 
showed that the differential proteins were mainly distributed 



in response to stimulus (50), metabolic process (32), locali- 
zation (7), and signalling (5), indicating that these processes 
play a leading role in fruit ripening (Fig. 8A). Catalytic activ- 
ity (48), binding (16), structural molecule activity (10), and 
transporter activity (5) were the the most abundant molecular 
function categories, which implied that catalysis and transpor- 
tation of substrates are vital for the fruit ripening (Fig. 8B). 
For cellular components, the organelle-related component 
was the largest group of proteins (83) (Fig. 8C). These results 
showed that the distributions of the differential proteins in 
functional categorization were consistent with the DEGs at 
the transcription level, which were also mainly distributed in 
categories of metabolic process, response to stimulus, orga- 
nelle, and catalytic activity (Fig. 4). 

As shown in Supplementary Fig. SI available at JXB 
online, there were 14, 5, and 8 differential proteins that were 
involved in 'post-translational modification, protein turno- 
ver, chaperones' of the Clusters of Orthologous Groups 
(COG) proteins, at 170, 190, and 210 DAF, respectively, and 
there were more downregulated proteins than upregulated 
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Fig. 5. qRT-PCR validation of differential expression. (A, B) Transcript levels of 1 9 genes, of which 1 1 involve plant hormone synthesis and signal 
transduction (A), and eight are most likely associated with fruit development and ripening (B), in WT (solid lines) and MT (broken lines) and the 
corresponding expression data of RNA-seq. WT-R: the data for RNA-seq in WT; MT-R: the data for RNA-seq in MT. The y-axis shows the relative gene 
expression levels analysed by qRT-PCR. The bars represent SE (n=3). (C) Comparison between the gene expression ratios obtained from RNA-seq data 
and qRT-PCR. The RNA-seq log 2 value of the expression ratio (/-axis) has been plotted against the developmental stages (x-axis). 
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Fig. 6. Transcript levels of genes in the ABA biosynthesis pathway (A), ABA signal transduction pathway (B,) and ethylene biosynthesis and signal 
transduction pathways (C). Bars represent the SE (n=3). An asterisk represents a statistically significant difference (P<0.05) analysed using Student's f-test. 



proteins (16 downregulated, five upregulated), includ- 
ing seven HSPs (Cs7gl2130.1, Cs9gl9220.1, Cs8gl8020.1, 
Cs8gl9490.1, Cs2g24040.1, Cs8gl8360.1, Cslgl2560.1), 
two glutathione S-transferase (orangel.lt03730.1, 



Cs6g03820.1), one zinc finger protein (orangel.lt02815.1), 
one peroxiredoxin (Cs6gl 3880.1), two cysteine protease 
(Cs4g07410.1, Cs3g23180.1), five proteasome regulatory sub- 
unit (Cslg08770.1, Cs3gl4360.1, Cslg25690.1, Cs7g07630.1, 
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Cs8g 16970.2) and three other proteins (Supplementary Table 
S10 available at JXB online). Next, there were 12, 12, 12, 
and 9 differential proteins between MT and WT involved in 
'general function prediction only', 'carbohydrate transport 
and metabolism', 'energy production and conversion', and 
'translation, ribosomal structure and biogenesis', respectively 
(Supplementary Fig. S 1 and Table S 1 0 available at JXB online ) . 
Notably, some proteins have been reported to be related 
to ripening, such as pectin methylesterase (Cslgl6560.1, 
Cslgl6550.1), citrate synthase (orangel.lt01588.1), and 
malic enzyme (Cs4gl9200.1, Cs4gl5270.1), which were all 
identified as differential proteins in MT compared with WT. 

As the transcript data were obtained from exactly the same 
powdered fruit pulp samples, we determined the number 
of identified proteins for which corresponding transcripts 
were represented in the RNA-seq data. The distribution of 
the corresponding mRNA:protein ratios is shown by a scat- 
terplot analysis of the log 2 -transformed ratios (Fig. 9). Of 
the 130 identified differential proteins, 54 had correspond- 
ing transcripts in the RNA-seq data. As shown in Fig. 9, 
almost all of the mRNA:protein ratios were concentrated 
at the centre of the plot (quadrant e), where protein and 



mRNA levels did not vary above 1.5- and 2-fold, respectively. 
Off centre, a total of nine mRNAiprotein ratios across all 
three stages were found in which both the mRNA and pro- 
tein levels exceeded this level of variation (Table 4, Fig. 9). 
The 170 DAF stage was characterized by relatively more 
mRNA:protein ratios falling into the quadrants a, c, g, and 
i compared with the 190 and 210 DAF stages, where the 
mRNA:protein ratios were substantially different. Moreover, 
five downregulated proteins, which reflected significant 
downregulation at the transcript level, were detected, namely, 
HSP18.1A (Cs8gl8020.1), ferritin (Cs6g09150.1), early nod- 
ulin-like protein 1 (orangel.lt03414.1), 18.2kDa class I HSP 
(Cs2g24040.1), and xyloglucan endotransglucosylase/hydro- 
lase protein 10-like (Cs8g03550.1). However, contrasting 
levels were observed for one gene, which exhibited upregu- 
lated expression of 40S ribosomal protein S9 (Cs7g04020.1) 
but downregulation of its protein abundance. At 190 DAF, 
we observed that mRNA:protein ratios fell mainly in quad- 
rants b, e, and h. Therefore, only a few mRNA:protein ratios 
reflected significant changes at both the transcript and pro- 
tein levels. Only two proteins displayed upregulation at 
both levels: polygalacturonase precursor (Cs9g01630.1) and 
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and WT. The proteins were categorized based on GO annotation, and the 
lecular functions (B), and cellular components (C). (This figure is available in 



remorin (orangel.lt00340.1). The 210 DAF stage was asso- 
ciated with the fewest mRNA:protein ratios outside of the 
central quadrant, and only one protein fell in the g quadrant: 
Cu/Zn superoxide dismutase (Cs3gl2080.1). 

As shown in Fig. 9, a total of 59 mRNA:protein ratios 
across all three stages were found to fall in quadrants b, d, 
h, and f, where the mRNA:protein ratios reflected signifi- 
cant changes at one level, the transcript or protein level, and 
showed a poor correlation between the transcripts and pro- 
teins. In summary, based on these results, a substantial degree 
of post-transcriptional regulatory activity during citrus fruit 
maturation is proposed, which has not been described in 
previous studies that focused solely on transcript or protein 
analysis. 

Discussion 

Fruit ripening is a highly coordinated, genetically pro- 
grammed, and irreversible phenomenon that involves a series 



of physiological, biochemical, and organoleptic changes that 
lead to the development of a soft, edible ripe fruit with desir- 
able attributes (Prasanna et al, 2007). In the present study, 
RNA-seq and iTRAQ technologies were used to investigate 
the differences in the transcriptome and proteome between 
the late-ripening bud MT and its WT. Hundreds of genes 
that were differently regulated during the three stages of cit- 
rus fruit ripening were identified by transcriptomic profiling. 
A total of 130 proteins identified by iTRAQ showed varia- 
tions in abundance during ripening in MT when compared 
with WT. The stage 170 DAF was associated with the largest 
number of differently expressed proteins with respect to WT. 
Of the 130 identified varying proteins, 54 corresponded to 
gene sequences that were obtained by RNA-seq, which ena- 
bled a comparison of ripening-related differences in specific 
transcripts or cognate proteins. In the present study, 68.5% 
(37 out of 54) of the protein and transcript pairs decreased 
or increased in parallel during ripening. This correlation was 
especially apparent at the early ripening stage (170 DAF). 



1666 I Wuefa/. 



170 DAF 190 DAF 210 DAF 




Transcripts Transcripts Transcripts 



Fig. 9. Comparison of changes in mRNA and cognate protein abundance. The relative change in abundance (MTVWT) is shown on a log 2 scale from 
samples from 170, 190, and 210 DAF fruits. Each colour denotes an mRNA:protein ratio: brown, 40S ribosomal protein S9 (Cs7g04020.1); orange, 
HSP1 8.1 A (Cs8g1 8020.1); blue, ferritin (Cs6g09150.1); yellow, early nodulin-like protein 1 (orange1.1t03414.1); red, 1 8.2 kDa class I heat-shock protein 
(Cs2g24040.1); light blue, xyloglucan endotransglucosylase/hydrolase protein 10-like (Cs8g03550.1); green, remorin (orangel .1t00340.1); deep red, 
polygalacturonase precursor (Cs9g01 630.1); and pink, Cu/Zn superoxide dismutase (Cs3g1 2080.1). (This figure is available in colour at JXB online.) 



A large number of genes showed consistency between the 
transcript and protein levels, whereas those genes that exhib- 
ited inconsistency between these levels suggest that post- 
transcriptional regulation plays an important role in the 
regulation of fruit ripening. Both the transcriptomic and 
proteomic data are important in deciphering the molecular 
processes involved in fruit ripening. The integrative transcrip- 
tomic and proteomic data not only highlighted a set of genes/ 
proteins that were possibly involved in the late-ripening trait 
formation of MT but also revealed molecular characteriza- 
tions associated with fruit development and ripening in sweet 
oranges. 

Identification of potential regulators and metabolism 
pathways involved in fruit ripening 

Some of the regulated genes and proteins in the most impor- 
tant pathways are shown in Fig. 10. Among these, the ABA 
pathway may play a central role in orange fruit development 
and ripening, and this regulation may function in combi- 
nation with other hormones, including ethylene. We also 
observed upstream and/or downstream biological processes, 
such as sugar metabolism and cell-wall biosynthesis, that 
were significantly changed during fruit ripening. These data 
depict a detailed picture of the regulation network involved 
in orange fruit ripening. 

The phytohormone ABA is a major endogenous factor 
and is the primary signal that regulates fruit development, 
maturation, and senescence (Gosti et al, 1999; Romero et al, 
2012; Yun et al, 2012; Jia et al, 2013b). In non-climacteric 
fruit, ripening is thought to be regulated by ABA, with eth- 
ylene also playing role (Koyama et al, 2010; Ji et al, 2012; 
Paul et al, 2012). The ABA level in plants is regulated by 
the pathways of ABA biosynthesis and degradation. NCEDs 
and A AOs are the key rate-limiting steps in ABA biosyn- 
thesis; ABA S'-hydroxylase is a key gene for the catabolism 
of ABA (Seo et al, 2000; Rodrigo et al, 2006; Zhang et al, 
2009). In this study, NCEDs, AAOs and ABA 8'-hydroxylase 



genes were found to be differentially expressed between MT 
and WT (Figs 6 and 9, Table 2). In addition, these genes were 
also identified to be differentially expressed in different rip- 
ening stages during fruit ripening in MT and/or WT (Figs 6 
and 9, Table 3). Furthermore, the curve of ABA of MT was 
different from that of WT (Fig. 1C). This, we hypothesized 
that ABA is a key factor for the formation of late ripening 
in MT, and that the genes involved in ABA biosynthesis and 
signal transduction pathways might be the mutation source. 
Two full-length cDNAs (CsNCEDl and CsNCEDl) encod- 
ing NCEDs were isolated and characterized from the epicarp 
of the orange fruits. CsNCEDl is likely to play a primary role 
in the biosynthesis of ABA, and CsNCEDl appears to play 
a subsidiary role that is restricted to chromoplast-containing 
tissue (Rodrigo et al, 2006). By suppressing the expression 
of SINCED1 using an RNAi method in tomato, ABA accu- 
mulation and SINCED1 transcript levels were downregulated 
to between 20 and 50% of the levels measured in the con- 
trol fruit, which led to downregulation in transcription of the 
genes that encode major cell- wall catabolic enzymes, such as 
PG, PME, XET, and others (Sun et al, 2012). Accordingly, 
in our study, the NCED1 gene was also identified as having 
differential expression during fruit ripening between MT and 
WT (Fig. 6). Thus, NCED1 might be a significant regulator 
for orange fruit ripening. The initial response to ABA implies 
the ABA-dependent PYR/PYL-mediated inactivation of 
PP2Cs, which allows for the release of SnRK2s and, hence, the 
phosphorylation of ABA-dependent transcription factors or 
other proteins. This ABA signal is later attenuated by upreg- 
ulation of the PP2Cs and the downregulation of the PYR/ 
PYL genes in an ABA-dependent manner. Thus, resetting of 
the ABA transduction pathway offers a dynamic mechanism 
for modulating the ABA response (Nishimura et al, 2010). 
Interestingly, in the present study, upregulation of the PP2Cs 
genes (such as CsABIl, CsABI2, and CsAHGl, among oth- 
ers) was concomitant with downregulation of the PYL4 
gene in both WT and MT (Figs 5 and 6, Table 3). Therefore, 
these results suggest that a transcriptional negative-feedback 
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regulatory mechanism might modulate the ABA responses 
during citrus fruit ripening, which is consistent with a previ- 
ous study (Romero et al, 2012). At the 190 DAF stage, we 
found that remorin protein was upregulation in MT at tran- 
script and protein levels (Table 4, Fig. 10). In a previous study 
(Lin et al, 2003), remorin homologues could be upregulated 
by exogenous ABA, which revealed that remorin protein 
was involved in the ABA signal transduction pathway. In 
addition, experiments in vivo and in vitro have proved that 
oligonucleotide galacturonic acid can significantly improve 
the phosphorylation level of remorin protein, and oligo- 
nucleotide galacturonic acid could influence the triggering 
of some physiological processes that depend on hormones 
(Bellincampi et al, 1996). 

Ethylene is also involved in the ripening of orange fruits, which 
are non-climacteric, and has long been known to play a major 
role in the ripening process of climacteric fruits (Theologis, 1992; 
Alexander and Grierson, 2002; Klee, 2002). Increasing evidence 
has indicated that ethylene is also involved in regulation of the 
ripening of non-climacteric fruits; even the small amount of 
ethylene that is produced by a fruit might be sufficient to trig- 
ger ripening-related physiological responses (Perkins- Veazie 
et al, 1996; Chervin et al, 2004; Trainotti et al, 2005; Iannetta 
et al. , 2006). Experiments with 1-methylcyclopropene (an action 
inhibitor of ethylene) have indicated that ethylene is required 
for the onset of accumulation of anthocyanins, fruit swelling, 
and the decrease in acidity that is associated with the ripening 
of grape berries (Chervin et al. , 2004). In the present study, two 
genes were identified as ACO; these were downregulated in MT 
(Fig. 10, Table 2). ACO is a key enzyme for ethylene biosynthe- 
sis. Two ETR, one EBF4, and one ERF1B genes were identified 



as being involved in the ethylene signal transduction pathway; 
these genes were also downregulated in MT (Fig. 10, Table 2). 
Ethylene is perceived by the receptors ETR1 and related proteins 
(Chang et al, 1993; Hua and Meyerowitz, 1998). The ethylene 
signal is transduced to EIN3 through CTR1 (Kieber et al. , 1993) 
and EIN2 (Alonso et al, 1999). Studies have also revealed that, 
in the absence of ethylene, the EIN3 protein is quickly degraded 
through a ubiquitin-proteasome pathway that is mediated by 
two F-box proteins (EBF1 and EBF2), whereas the EIN3 protein 
is stabilized by ethylene itself (Guo and Ecker, 2003; Potuschak 
et al, 2003; Gagne et al, 2004). Citrus (and other non-climac- 
teric fruits) lack a ripening-associated autocatalytic rise in eth- 
ylene production. However, Goldschmidt (1997) pointed out 
the potential significance of even the lowest levels of ethylene 
present in non-climacteric fruits in regulation of the ripening 
process. Citrus reveal ripening-related symptoms in response to 
exogenous ethylene, and induction of the ripening of citrus fruit 
peel by ethylene was also found to be opposed by gibberellins 
and cytokinins, possibly due to the reduction in the tissue's sensi- 
tivity to ethylene (Goldschmidt, 1997). Recent work on apple by 
Johnston et al. (2009) has revealed that early- and late-ripening 
events differ not only in their dependence on ethylene but also 
in their sensitivity to ethylene. In the present study, we observed 
that some genes that are related to ethylene biosynthesis and 
signal transduction were differentially expressed between MT 
and WT, and were also differentially expressed in the different 
fruit-ripening stages of WT and MT (Fig. 10, Tables 2 and 3). 
These findings thereby necessitate further work not only on the 
role of the development/developmental signals but also on the 
mechanism of alteration in the sensitivity of the tissue/s with the 
development of fruit. 




Sugar and acid metabolism 



Fig. 10. Summary of some of the biological pathways involved in citrus fruit ripening. Red boxes indicate genes/proteins that were upregulated in 
MT compared with WT, brown boxes indicate genes/proteins that were upregulated during fruit ripening in WT, green boxes indicate genes/proteins 
that were downregulated in MT compared with WT, grey boxes indicate genes/proteins that were downregulated during fruit ripening in WT, and blue 
boxes represent genes/proteins that were shown to have mixed expression patterns (both up- and downregulation) in MT compared with WT. ACO, 
1-aminocyclopropane-1-carboxylate oxidase; ETR, ethylene receptor; EBF, EIN3-binding F-box protein; ERF1 , ethylene-responsive transcription factor; 
NCED, 9-cis-epoxycarotenoid dioxygenase; AAO, abscisic-aldehyde oxidase; PYR/PYL, abscisic acid receptor; PP2C, protein phosphatase 2C; 
SnRK2, serine/threonine-protein kinase SRK2; XET, xyloglucan endotransglycosylase; PG, polygalacturonase; PE, pectinesterase; GAUT, alpha-1 ,4- 
galacturonosyltransferase; SPS, sucrose-phosphate synthase; SS, sucrose synthase; CS, citrate synthase; and IDH, isocitrate dehydrogenase. (This 
figure is available in colour at JXB online.) 
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Sugars, as crucial components of fruit quality, are also 
involved in signal transduction. The effect of sugar on many 
developmental processes is concentration dependent (Gibson, 
2005). In the present study, the sugar content of MT was lower 
than that of WT (Fig. 2D-F); however, the citric acid content 
in MT was higher than that in WT (Fig. 2B). Furthermore, 
some genes/proteins that are involved in sucrose, citrate, and 
fructose metabolism were downregulated in MT (Fig. 10, 
Table 2). Interestingly, in the present study, we found that 
changes in the sucrose content were closely correlated with 
phase transitions in the developmental process that spans the 
mature green phase to full ripening; the content of sucrose 
remained low and relatively stable in the early stages of the 
mature green phase and increased sharply when the fruits 
began to enter the break phase (Fig. 2F). In the strawberry, 
sugars, especially sucrose, can promote the ripening of fruit 
(Jia et al, 2011). In the present study, a citrus sucrose trans- 
porter 1 gene was identified in MT and WT; this gene was sig- 
nificantly changed during fruit ripening (Table 3). The RNAi 
of the sucrose transporter 1 (SUT1) in the strawberry was 
able to significantly reduce the sucrose content and delayed 
fruit ripening (Jia et al, 2013a); therefore, sucrose deserves 
further investigation as a potentially significant regulator for 
citrus fruit ripening. 

Cell-wall metabolism and other gene/protein changes 
during fruit ripening 

Cell-wall metabolism is one of the major fruit ripening- 
related changes; this type of metabolism involves the disman- 
tling of multiple polysaccharide networks by diverse families 
of cell-wall-modifying proteins, including enzymes for pectin 
and cellulose catabolism. The degradation of pectin and cel- 
lulose depends on ethylene during the softening of climac- 
teric fruits (Ergun et al. , 2005; Nishiyama et al , 2007). Recent 
work has revealed that ABA affects cell-wall catabolism dur- 
ing fruit ripening via downregulation of the expression of 
major catabolic genes {PG, PE, and XET) (Sun et al, 2012). 
A subset of cell-wall-modifying genes and proteins was found 
to be downregulated in MT and WT during fruit ripening, 
including PEs, PGs, and GAUTs (Fig. 10, Table 3). GAUT1 is 
the core of synthesis of homogalacturonan, the most abun- 
dant pectic polysaccharide and the core of the plant cell- 
wall pectin biosynthetic homogalacturonamgalacturonosylt 
ransferase complex (Atmodjo et al, 2011). Comparing MT 
with WT, the genes/proteins that were identified as PGs were 
upregulated in MT and (3-glucosidase genes, and XETs dis- 
played mixed expression patterns (Fig. 10, Table 2). Changes 
in the activity of these cell-wall-related genes are known to 
result in the abnormal development of juice sac granulation 
(Sharma and Saxena, 2004), and modifications in the cell- 
wall structure or in the components of the membranes of the 
segments and juice sacs during fruit development and ripen- 
ing clearly influence the formation of the fruit pulp melting 
characteristic (Waldron et al, 2003). 

A striking feature of the proteomic data was that seven 
proteins (Supplementary Table S10 available at JXB online) 
were identified as HSPs, and 1 1 differential expression HSP 



transcripts were also identified in MT during fruit ripening 
(Table 3). The possible important roles of HSPs in fruit devel- 
opment and ripening have been reported recently in tomato 
(Neta-Sharir et al , 2005; Faurobert et al. , 2007), apple (Wang 
et al, 2009), apricot (Grimplet et al, 2005), and citrus (Pan 
et al, 2012). In the present study, seven HSPs, including two 
HSP90s exhibiting increased the abundance in MT and five 
small HSPs exhibiting decreased the abundance in MT were 
identified (Supplementary Table S9 available at JXB online). 
HSPs, which function as chaperone molecules, are important 
in post-translational modification, and the seven HSPs iden- 
tified in this study indicate possible roles in fruit ripening. 

In the present study, in addition to ABA and ethylene, 
sucrose appeared to function as another key signal in the regu- 
lation of citrus fruit ripening. Fruit ripening is the result of an 
orderly alteration of a series of physiological and biochemical 
events, such as sugar and acid metabolism, cell-wall metabo- 
lism, and pigment and flavour metabolism, and each of these 
metabolic systems or processes is involved in the regulation of 
a number of genes. It is not known whether each metabolic 
system is independently regulated by a unique signalling path- 
way or by a common signalling network. It is unlikely that 
many metabolic systems across a number of developmental 
stages would be controlled by a single signal. Thus, it is not 
surprising that the development and ripening of citrus fruit 
may be controlled by multiple signals. However, if one signal 
pathway changes, then this change could affect the entire pro- 
cess of fruit development and ripening. Previous studies on 
strawberry showed that exogenous ABA promoted strawberry 
ripening while repressing ABA levels by the administration 
of an ABA biosynthesis inhibitor, which resulted in a delay 
in ripening and a decrease in the sugar content; in addition, 
treatment with glucose and especially sucrose increased the 
levels of ABA and promoted fruit ripening (Jia et al, 2011, 
2013a). These data suggest that the two signalling pathways 
interact cooperatively to regulate fruit ripening. In the pre- 
sent study, the ABA pathway, ethylene pathway, and sucrose 
pathway were all different between MT and WT, and we have 
not proven which is the most important factor for forming the 
late-ripening trait of MT. This issue will require further study. 



Conclusions 

A multiple-level analysis of the gene/protein expression 
changes in MT, compared with WT, was conducted in our 
study. Combined with transcriptomic and proteomic data, 
the results revealed several key candidate regulators that 
may play high-priority roles in citrus fruit ripening. Based on 
an integrated analysis of the transcriptome and proteome, 
we identified biological processes that appear to be of high 
importance in the ripening process, such as the signal trans- 
duction of hormone, the regulation of hormone and sugar 
levels, and cell-wall metabolism. Thus, the integrated analysis 
enabled us to generate additional information for a compre- 
hensive understanding of biological events that are relevant 
to the regulation of citrus fruit ripening. These results have 
revealed multiple ripening-associated events during citrus 
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ripening, providing new insights into the molecular mecha- 
nism of citrus ripening regulatory networks. 
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